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ABSTRACT: Previously identified QTL regions on BTAl and BTA5 were investigated to validate the QTL regions and to identify 
candidate genes for growth and carcass traits in commercial cattle populations from the USA and Korea. Initially, a total of 8 
polymorphic microsatellite (MS) markers in the BTAl and 5 QTL regions were used for Chi-square tests to compare the frequencies of 
individual alleles between high and low phenotypic groups for the US (Michigan Cattleman's Association/Michigan State University; 
MCA/MSU) cattle. For a subsequent study, 24 candidate genes containing missense mutations and located within the QTL regions based 
on bovine genome sequence data were analyzed for genotj^ping in the two commercial cattle populations. Re-sequencing analyses 
confirmed 18 public missense SNPs and identified 9 new SNPs. Seventeen of these SNPs were used for genotyping of the MCA/MSU 
cattle (n = 98) and Korean native cattle (n = 323). On BTAl, UPKIB, HRG, and MAGEFl polymorphisms residing between BM1312 
and BMS4048 were significantly associated with growth and carcass ttaits in one or both of the MCA/MSU and Korean populations. On 
BTA5, ABCD2, IL22 and SNRPF polymorphisms residing between BL4 and BR2936 were associated with marbling and backfat ttaits in 
one or both of the MCA/MSU and Korean cattle populations. These results suggested that BTA 1 and 5 QTL regions may be segregating 
in both Korean Hanwoo and USA commercial cattle populations and DNA markers tested in this study may contribute to the 
identification of positional candidate genes for marker-assisted selection programs. (Key Words: Quantitative Trait Loci, Microsatellite 
Markers, SNPs, Carcass Traits) 



INTRODUCTION 

Quantitative trait loci (QTL) analyses have been 
performed for economically important quantitative traits in 
beef cattle, since such genomic technologies can provide 
genetic dissection of polygenic traits and help to understand 
the molecular mechanisms of the identified genes or genetic 
loci (Womack, 2006). However, effectiveness of these QTL 
must be validated in commercial populations before 
implemented into marker-assisted selection. Validation of 
such QTL effects in commercial populations includes i) the 
size of the effects should be sufficiently large to be 
detectable with a reasonable size of sample populations 
even under non-genetic or environmental influences; ii) the 
pattern of the effects should be predictable from the 
selection even in different genetic backgrounds (or lines), 
iii) the selection of the favorable allele for the trait of 
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interest should not have a deleterious effect on other traits. 

In an experimental population derived from a Bos 
indicusxBos taurus cross, Kim et al. (2003) reported a QTL 
in the interstitial region of bovine chromosome 1 (BTAl) 
between BM1312 (61.85 cM) and SST (80.38 cM) affecting 
hot carcass weight, slaughter weight, and yearling weight. 
Later Morsci et al. (2006) reported the presence of QTL 
affecting marbling, ribeye muscle area and fat thickness in 
the vicinity of SST on BTAl in Angus, but excluded SST as 
underlying the ribeye muscle area QTL. Also QTL for fat 
yield have been mapped to the intermediate region of BTAl 
between BMS945 (51 cM) and n.STS083 (67 cM) in a 
crossbred cattle population (Casas et al., 2003; Casas et al., 
2004). 

On BTA5, three chromosomal regions (0 to 30 cM, 55 
to 70 cM and 70 to 80 cM) were identified as having 
significant associations with growth traits (Li et al., 2002), 
and two backfat QTL were mapped to the intermediate 
region of BTA5 (RM500 to BR2936 and BMS490 to 
ETHIO) in a commercial Bos taurus line (Li et al., 2004 j. A 
QTL region of BTA5 influencing retail fat contents, fat 
yield and fat thickness in half-sib or full-sib progeny from 
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different cattle breeds (Bos taurus) was also reported by 
Casas et al. (2000). According to these studies, both BTA 1 
and 5 QTL were associated with several beef growth and 
carcass traits, including birth weight, yearUng weight, 
slaughter weight, hot carcass weight, adjusted backfat 
thickness, average daily gain, marbling, texture and 
tenderness. 

Therefore, the objectives of this study were i) to vaUdate 
the segregation of the QTL on BTAl and 5 in different 
commercial cattle populations without pedigree information 
and ii) to further characterize these QTL regions using 
microsatellite markers and single nucleotide polymorphisms 
(SNPs) from candidate genes for growth and beef quality 
traits. 

MATERIALS AND METHODS 

Animals and phenotypic data 

Commercial cattle from the USA used for this study 
were a part of the Michigan Cattlemen's Association 
(MCA)/Michigan State University (MSU) Farm to Fab 
Program that was a means for Michigan beef producers to 
obtain feedlot performance and carcass merit information 
for a portion of their calves. For the present study, steers (n 
= 100) representing eight breeds enrolled in the Farm-to- 
Fab Program in the fall of 1997 were used (Angus, Polled 
Hereford, Maine Anjou, MARC 11, Salers, Shorthorn, 
Simmental and Tarantaise). End point was established as 
0.89 cm of backfat, measured by ultrasound, or 635.2 kg 
live weight, whichever occurred first. Phenotypic traits 
recorded included birth weight (BW), average daily gain 
(ADG), back fat thickness (BF), slaughter weight (SW), 
marbling score (MAR), Warner-Bratzler shear force 
(WBSF), and ribeye area (REA). DNA was isolated from 
white blood cells using standard procedures. 

Korean cattle phenotypic data for carcass and beef 



quality traits were obtained from 323 Hanwoo steers (30 
months of age) slaughtered at Pyungnong slaughterhouses, 
Korea between October and December 2007. The measured 
phenotypes according to the Korean Beef Marbling 
Standard were carcass weight (CW), longissimus muscle 
area (LMA), backfat thickness (BFT), and marbling score 
(MS), but their feeding program, concentration and forage 
intake and fattening period were not available. Genomic 
DNAs were extracted using a QIAprep® Spin Miniprep Kit 
according to the manufacturer's protocol (Qiagen, USA). A 
description of phenotypes measured for this study is shown 
in Table L 

Genotyping of microsatellite markers 

Bovine chromosomes 1 and 5 were chosen for their 
association with carcass and growth traits (Li et al., 2002; 
Kim et al., 2003). Four microsatellite (MS) markers 
spanning QTL regions on each chromosome were 
genotyped in the MCA/MSU cattle. Selection of the MS 
markers was based on the relative position and fragment 
size from the genomics database at the USDA, ARS, US 
Meat Animal Research Center (MARC) (http://www.marc. 
usda.gov/). Primer sequences and amplification conditions 
were also obtained from the MARC database. The BTAl 
MS markers were BM1312, BMS4001, BM6506 and 
BMS4048. The BTA5 MS markers were BL4, BR2936, 
ETHIO and BMS1248. PCR products were separated by 
polyacrylamide gel electrophoresis with an ABI 373 
Automated Laser Fluorescent DNA sequencer (PE Applied 
Biosystems, Foster City, CA, USA), and aUeUc sizes were 
determined with GeneScan® software. 

Single nucleotide polymorphism detection 

A total of 25 pairs of primers were designed for 
amplification and sequencing of missense SNPs in 24 
candidate genes located in the growth and carcass trait QTL 



Table 1. Means and standard deviations for growth, carcass and beef quality traits measured on MCA/MSU and Korean cattle 
populations 

Phenotypes N Mean StdDev Median 

Michigan cattle 

Average daily gain (ADG, kg) 100 

Back fat thickness (BF, cm) 97 

Birth weight (BW, kg) 91 

MarbUng (MAR) 97 

Warner-Bratzler Shear force (WBSF, kg) 97 

Slaughter weight (SW, kg) 1 00 

Ribeye area (REA, cm^) 97 
Korean cattle 

Carcass weight (CW, kg) 319 

Longissimus muscle area (LMA, cm^) 320 

Backfat thickness (BF, cm) 321 

MarbUng score (MAR) 1-7 320 



1.37 
0.97 
40.86 
530.52 
3.39 
579.46 
34.321 

346 
79.3 
5.04. 
2.52 



0.19 

0.2 

8.07 

71.57 
0.78 

51.29 
3.63 

38 
8.31 
2.93 
1.30 



1.37 
0.89 
39.92 
520.00 
3.30 
579.92 
34.79 

344 

79 
4 
2 
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regions on BATl and BTA5. These genes included 11 genes 
on BTAl spanning the region between the BM1312 and 
BMS4048 markers, and 13 genes on BTA5 spanning the 
region between the BL4 and BMS1248 markers. Primer 



information for PCR and sequencing is listed in Table 2. 
Polymerase chain reactions were performed in 10 yl 
volumes containing 12 ng of genomic DNA, 10 pmol of 
each primer, 200 \iM of each dNTP, 2.5 units of HS-Taq 



Table 2. Identification of SNPs in MCA/MSU and Korean cattle 



Chr. 



Gene 
symbol 



GenBankJD 



Primer sequences 
(5'^3') 



Product 
size (bp) 



NCBl' 



Site of SNPs in this study 



UPKIB rs43652277 TTTGGTTGCCTGGCAGTGAG 321 

ATGCCTACGATGCCCAGGAC 
PARP14 rs43242284 AATCCTGGATGAATGCTTTT 455 

ATTAACTGGTCCTTGGGAGA 
MYLK rs43244443 GAGGAGAGGCGGGGGTGCTG 435 

GGAGCAGGGGCCCAGGAAGA 
ITGB5 rs41257187 TCACGGTCCTCAGGGAGCCA 443 

CAGCCAGCCCTTCCTGCACA 
SST rs\lKlQ991 CCTCATTTCATCCTGCTCAG 307 

ATTTTGTGGCTGTAGCAACA 
KNG rs43705174 TCCTTCTCGCAGAAATGTG 449 

TGAACTGATCGGAAAGGTG 
HRG(l) rs43245364 AGCCCTCTTCACTTTGTCCA 336 

GCTCTGCCACTGATTGACTG 
HRG(2) rs43245367 CAACCACCCTGGGCAAGATA 344 

GGTAGGGCACCTGCAATGAC 
MAGEFl rs 17870725 CGGCCCTTGCGATGATCTCA 292 

TGACCGCCTCGCAGGAGATG 
B3GNT5 rs42820807 TGTTTGTTAGTGGCAGGAGA 419 

ATCAGTGGGTCAGAAGGTGT 
ABCC5 rs43247946 GGAGTGGCAGGAGCATGAT 426 

GTGCATCCCCACAGAAACC 
AP2M1 rs43246948 GGGCAGGGCTGTGGAACTTC 329 

GCAAAGGGCAGGCAGCAAAG 
5 ABCD2 rs42920919 AGGTGGCAGGGATGGCAATC 303 

CGGCTGCGGAGAACAGAGAA 
5 IL22 rs29()03650 CCAGAGTGAAGGTGTGGTTG 327 

ACGAACTTGGTAGGATTGGC 
5 DGKA rs43431357 TCCCTGAGAATTTGTGGAG 438 

CCAGTTTGGGAGTGATACA 
5 SNRPF rs29004136 TCCCTGAGAATTTGTGGAG 431 

CCAGTTTGGGAGTGATACA 

5 AMDHDl rs29002687 GGGCGCGCCTCAGACAACCA 384 

GGCCGCGCGCTACCAGGAAG 
5 RDH16 rs41590097 CCCAGGGAGCTGCCAACCAG 371 

GGTGGTCCTCGTGGGCCTGT 
5 PPT rs43087855 ATTCCAGCATATTCCTGACA 451 

CCTGCATAGGAGATGGACTT 
5 SPIC rs43440868 TTCCAAGAAATAAGGTGCAG 382 

CACTCAACAAATCCTTCCTG 
5 DRAM rs42500325 TCGGGCTTGCAAGGACCTGA 324 

GGGCGGGTGCATTTGTGATG 
5 HMOXl NM_001014912 TGAGGTCCGGCAGAGGGTCC 387 

CTGCCCCCTGCCCCAAGCTC 
5 TXN2 NM_174208 GGAGGGACCCGCTGTGCATC 365 

TGCTGCTTGGCCACCACCTT 
5 CACNG2 rs43436464 GGTCCTGCGGTTGGCTGTGT 489 

TCTGCCAATGCCGGAGACCC 
5 C1QTNF6 rs43438424 TTCCAGCCGCTGCTCTTTGA 464 

ACCAGGTCAGGGGCCCATTT 



25 IR 

25 IR 

25 lY 

27Y 

177Y 

lOlS 

25 IM 

251Y 
256S, 300K 

251S 

25 IW 

251K 

25 IM 
73R, 158Y 

251S 

18Y,21R, 113M, 
138M, 178M, 
256R, 258S 
256R, 257Y, 304R 

251Y 

251Y 

251Y 

251S 
Unknown 
Unknown 

25 IM 

25 IR 



146(C/G)'', 190 {Cn)\209{PJGy- 
251(A/G)',346(A/G)' 
251(A/G)'-' 



251(A/C)'-^ 
251(C/T)'-' 
256(C/G)', 300(G/T)'-' 



251(A/C)''' 
73(A/G)'•^ 158(C/T)' 



113(A/C)'', 138(A/C)'- 
178(A/C)', 258(C/G)'-- 



304(A/G)' 
251 (C/T)' 



324(C/T)^- 



135(A/G)^ 168(A/G)''\ 229{CrTf-\ 
251(A/G)''\ 282(A/T)''' 



' SNP previously reported in the NCBI database and position is based on the reference sequence (http://www.ncbi.nlm.nih.gov/projects/SNP/). 

^ New SNPs based on bovine SNPs at the reference sequence. 

' Seventeen SNPs that were genotyped in the MCA/MSU cattle and Korean cattle for association analyses. 
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DNA polymerase (Solgent , Korea), and reaction buffer 
with 1.5 mM MgCla. Thermocycling was performed in a 
PTC -200 thermocycler (MJ Research, Watertown, MA, 
USA) with a 10 min initial denaturation at 94°C, 35 cycles 
of 94°C for 30 s, 45 to 62°C for 40 s, 72°C for 30 s, and a 
final extension at 72°C for 10 min. The PGR products were 
identified by 2% agarose gel electrophoresis at 100 mV for 
20 min and sequenced by Solgent (Taejeon, Korea) with an 
ABI PRISM 3730 DNA Analyzer (Applied Biosystems, 
Foster City, CA, USA) according to the manufacturer's 
instructions. The sequencing results were analyzed for 
polymorphism detection using the Sequencher program 
(Gene codes, version4.7, Ann Arbor, MI, USA). 

Genotyping by single-base extension (SBE) and 
electrophoresis 

For genotyping SNP sites, amplification and extension 
primers were designed for the single-base extension (SBE) 
method (Vreeland et al., 2002). Primer extension reactions 
were performed with the SNaPshot ddNTP Primer 
extension Kit (Apphed Biosystems). To clean up the primer 
extension reaction, one unit of SAP (shrimp alkaline 
phosphatase) was added to the reaction mixture, and the 
mixture was incubated at 37°C for 1 h, followed by 15 min 
at 72°C for enzyme inactivation. The DNA samples, 
containing extension products, and GeneScan 120 Liz size 
standard solution was added to Hi-Di formamide (Applied 
Biosystems) according to the recommendations of the 
manufacturer. The mixture was incubated at 95 °C for 5 min, 
follow by 5 min on ice, followed by electrophoresis using 



the ABI Prism 3100 Genetic Analyzer. The results were 
analyzed using the ABI Prism GeneScan and Genotyper 
software (Applied Biosystems). Probe information is 
available in Table 3. 

Statistical analysis 

Evaluation of the QTL regions on BTA 1 and 5 for both 
the MCA/MSU and Korean cattle populations were 
performed with several statistical methods that would 
account for genetic associations of marker alleles within the 
QTL regions. 

Significant differences of MS alleUc frequencies were 
determined fiom 2x2 contingency tables (high vs low 
phenotypic groups x with vs without particular allele) in 
terms of Chi-square. Other tests of 2x2 contingency such as 
continuity-adjusted chi-square and Fisher's exact test were 
also performed. To maximize the sample size for allelic 
frequency comparison, medium phenotypic values were 
used to divide the high and low phenotypic groups in the 
MCA/MSU cattle. 

GLM procedures of SAS (Version 9.01; SAS, Inst., Inc., 
Gary, NC, USA) were used to examine associations 
between SNP and phenotypes for both the MCA/MSU and 
Korean data sets. The general linear model used for the 
MCA/MSU data was: 

Y = n-HSBi-HG,-He 

where |i = overall experimental mean, SB, = fixed effect of 
the fth sire breed, G^ = fixed effect of genotype, and e = 



Table 3. List of probe information and allele frequencies of 17 polymorphisms on BTAl and BTA5 genotyped in MCA/MSU and 



Korean cattle 


Gene 


SNP_ID 


SNaPshot primer sequences 


'Allele r 


MSU allele 1 

frequency 


KOR allele 1 
frequency 


UPKIB 


209-R 


TTCAGTCAGGCTGATGGCC (F) 


A 


0.82 


0.71 


UPKIB 


146-S 


TGCTTACAGCCTATGCAG (R) 


G 


0.79 


0.33 


PARP14 


251-Y 


TGGCAGTGCAGATTGCAGC (R) 


C 


0.68 


0.44 


HRG(2) 


251-Y 


CTGGGACGACTGTGAGCTG (R) 


C 


0.83 


0.91 


HRG(l) 


251-M 


GCTCAGTATCCTCAAAGAAAT (F) 


A 


0.78 


0.86 


MAGEFl 


300-K 


AGCAGAAACTGCACCAACTC (F) 


G 


0.59 


0.81 


ABCD2 


251-M 


TCTGCATTCACTCCAGGTGA (F) 


A 


0.59 


0.74 


IL22 


73-R 


CGTTGAGCTTACAGTGCGAC (F) 


G 


0.83 


0.72 


SNRPF 


113-M 


GGATGAATGGGAAACGCGC (F) 


C 


0.52 


0.50 


SNRPF 


138-M 


ACCACCAGGCCCGGCAAGT (F) 


C 


0.87 


0.96 


SNRPF 


258-S 


GCAGGATACTCGCCGCCGA (F) 


C 


0.76 


0.98 


AMDHDl 


304-R 


GGAGTTGGACATCGGCATC (F) 


G 


0.71 


0.84 


DRAM 


324-Y 


CGTTCGTCTGCGGCGTCGT (R) 


T 


1 


0.94 


CIQTNF6 


168-R 


TCGCTGTAGACGGCATTATC (R) 


G 


0.93 


0.82 


CIQTNF6 


229-Y 


AAGAGCCGTGCCCAGACG (R) 


C 


0.65 


0.56 


CIQTNF6 


251-R 


CATCACGCTCTGGCTCTGC (R) 


A 


0.83 


0.56 


CIQTNF6 


282-W 


CGCGTACAGGATGACGGC (R) 


A 


0.85 


1 



MSU = MCA/MSU cattle population (n = 98); KOR = Korean cattle population (n = 323). 
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residual error. For the Korean cattle population, Yy = |i+G, 
+e, where G, = fixed elTect of the j'th Genotype, and e = 
residual error. Depending on the trait under consideration, 
slaughter weight or birth month were included in both 
models as covariates. 

RESULTS 



Korean cattle, and the DRAM (324 Y) and CIQTNF6 (168R) 
genes in the MCA/MSU cattle (Table 3). Allele frequencies 
of the 17 SNPs were similar between the MCA/MSU and 
Korean cattle populations, although frequencies for a few 
such as UPK1B-U6S, PARP14-251Y and CIQTNF6-25m 
were found to be different between the Korean and 
MCA/MSU cattle (Table 3). 



SNP identiflcation and genotyping 

By direct DNA sequencing of 25 amplicons using 
random DNA samples from the MCA/MSU (n = 12) and 
Korean (n = 8) cattle, a total of 27 SNPs were identified 
from 12 primer sets, which included 18 public SNPs and 9 
new SNPs (Table 2). Fifteen previously reported SNPs from 
13 candidate genes in the BTAl and BTA5 QTL regions 
were not detected to be segregating in either the MCA/MSU 
or Korean cattle (Table 2). Based on these sequencing 
results, SNPs that met the following criteria were selected 
for genotyping: i) location (non synonymous 
polymorphisms were preferred), ii) minor allele frequency 
(MAF freq. >0.05), and iii) linkage disequilibrium (LD, 
only one polymorphism was chosen if it was in complete 
LD (r^ = 1) with one or more other polymorphisms). Pair- 
wise sequence analyses from the MCA/MSU and Korean 
cattle found that some polymorphisms within BTAl or 
BTA5 were in complete LD ( |D' = l| and r^ = 1). 

As a result, six SNPs {UPK1B-146S and 209R, 
PARPl 4-251Y, HRG(l)-25m, HRG(2)-251Y and 
MAGEF1-300K) in the BTAl QTL region, and eleven 
SNPs (ABCD2-251M, 7L22-73R, SNRPF-IUM, 138M and 
258S, AMDHDl -304R, DRAM-324Y and CIQTNF6-168R, 
229 Y, 25 IR and 282W) in the BTA 5 QTL region were 
selected to genotype the MCA/MSU (n = 98) and Korean 
native cattle (n = 323) using the SNaPshot method. 

Most genotyped SNPs had a MAF higher than 0. 1 with 
the exception of the SNPs in the SNRPF (138M and 258S), 
DRAM (324Y), and CIQTNF6 (282W) genes detected in 



Phenotypic association of BTA 1 and 5 markers with 
growth and carcass traits for MCA/MSU and Korean 
cattle 

Potential association of BTA 1 and BTA 5 QTL regions 
with growth and carcass traits in the MCA/MSU cattle 
population was initially assessed using MS markers. 
Individual alleUc frequencies of MS markers were 
compared between high and low phenotypic groups in the 
MCA/MSU cattle (Table 4), and chi-square tests were 
performed to determine significance levels. The BTAl MS 
markers exhibiting significant allelic frequency differences 
were BM1312 for Warner-Bratzler shear force and ribeye 
area (REA), and BM6506 for birth weight. Significant 
differences in allelic frequency for BTA5 markers were 
found with BL4 for average daily gain and birth weight, 
BR2936 for average daily gain, ETHIO for marbling and 
slaughter weight, and BMS1248 for average daily gain and 
slaughter weight. These results suggested that QTL for 
these growth and carcass traits were segregating in the 
MCA/MSU cattle. 

Association analyses were performed between SNP 
genotypes and growth and carcass traits using phenotypic 
data for the MCA/MSU and Korean catde (Table 5). The 
association analysis revealed that the UPK1B-2Q9R, 
ffl?G(l)-251M, and HRG(2)-251Y polymorphisms on 
BTAl were significantly associated with ribeye area (REA), 
marbling (MAR), and backfat thickness (BF) in the 
MCA/MSU cattle. The 146S SNP genotypes of the UPKIB 
gene were significantly associated with REA and MAR. 



Table 4. Chi-square analyses of allelic frequency differences between high and low phenotypic groups 



Chr. 



Marker (allele) 



Traits 



High' 



MCA/MSU cattle 



Low 



BTAl 


BM1312(109) 


REA* 


0.39 


0.61 


4.52 






WBSF* 


0.41 


0.59 


4.02 




BM6506(201) 




0.69 


0.31 


8.58 


BTA5 


BL4(159) 


ADG*** 


0.35 


0.65 


19.22 




BL4(165) 


BW* 


0.67 


0.33 


4.12 




BR2936(162) 


ADG* 


0.58 


0.42 


4.51 




ETH10(220) 


SW* 


0.35 


0.65 


4.50 




ETH 10(222) 


MAR*** 


0.67 


0.33 


8.99 




BMS1248(142) 


ADG** 


0.71 


0.29 


6.95 






SW* 


0.71 


0.29 


9.15 



' Medium phenotypic values were used to divide the high and low phenotypic groups 

* Significant at the 5% level. ** Significant at the 1% level. *** Significant at the 0.5% level. 
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Table 5. Association of SNPs on BTAl and BTA5 witli growtii and carcass traits for MCA/MSU and Koran cattle 



Chr. 


Gene 


Phenotypic 
trait 


MSU/MCA cattle (n 


= 98) 


p-value 


Phenotypic 


Korean cattle (n = 


323) 




Genotypic least squares means (SE) 


trait 


Genotypic least squares means (SE) 


- p-value 


BTAl 


UPKIB 


SW 


AA: 585.63 


AG: 563.52 


GG: 562.52 


0.0873 










NS 




(209A/G) 




(6.27) 


(12.07) 


(18.81) 


















REA 


AA: 34.45 


AG: 34.87 


GG: 31.02 


0.0527 


















(0.47) 


(0.886) 


(1.38) 














BTAl 


UPKIB 


RA 


CC: 37.23 


CG: 37.79 


GG: 33.66 


0.0422 














(146C/G) 




(0.74) 


(0.52) 


(0.49) 


















MAR 


CC: 585.714 


CG: 549.090 


GG: 518.727 


0.0297 


MAR 


CC: 2.897 


CG: 2.627 


GG: 2.976 


0.0497 








(16.556) 


(14.9oj) 


(9.477) 






(0.121) 


(0.123) 


(0.147) 




BTAl 


HRG(2) 


BF 


AA: 0.68 


AG: 0.95 


GG: 1.04 


0.0074 


CW 


AA: 330.333 


AG: 348.000 


GG: 360.020 


0.0167 




(251A/G) 




(0.02) 


(0.05) 


(0.11) 






(14.401) 


(5.375) 


(4.371) 




BTAl 


HRG(l) 


BF 


AA: 0.67 


AC: 0.95 


CC: 0.99 


0.0338 














(251A/C) 




(0.11) 


(0.04) 


(0.03) 














BTAl 


MAGEFl 










NS 


BF 


GG: 5.280 


GT: 4.566 


TT: 5.985 


0.0243 




(300G/T) 














(0.180) 


(0.230) 


(0.752) 


















MAR 


GG: 2.684 


GT: 2.770 


TT: 3.532 


0.0326 


















(0.106) 


(0.124) 


(0.323) 




BTA5 


ABCD2 










NS 


MAR 


AA: 2.899 


AC: 2.587 


CC: 2.307 


0.0045 




(251A/C) 














(0.106) 


(0.127) 


(0.217) 




BTA5 


IL22 


BF 


CC: 0.93 


CT; 0.94 


TT: 1.19 


0.0021 










NS 




(73A/G) 




(0.03) 


(0.04) 


(0.07) 














BTA5 


SNRPF 


MAR 


AA: 554.545 


AC: 508.709 


CC: 536.428 


0.0728 


MAR 


AA: 2.543 


AC: 2.782 


CC: 2.946 


0.0612 




(113A/C) 




(15.437) 


(13.004) 


(13.683) 






(0.138) 


(0.107) 


(0.155) 




BTA5 


C1QTNF6 


MAR 


CC: 511.621 


CT: 546.571 


TT: 547.272 


0.0832 


BF 


CC: 5.396 


CT: 4.723 


TT: 5.194 


0.0946 




(229C/T) 




(11.541) 


(11.866) 


(21.167) 






(0.242) 


(0.206) 


(0.315) 




BTA5 


C1QTNF6 


WBSF 


CC: 3.339 


CT: 3.822 




0.0969 










NS 




(168A/G) 




(0.094) 


(0.271) 

















Four polymorphisms {UPK1B-IA6S, MAGEF-3QQYi, 
ABCD2-25nA, and SNRPF-IUM) were significantly 
associated with MAR in the Korean cattle. Interestingly, the 
UPK1B-U(& and SNRPF-IUM polymorphisms exhibited 
significant association with MAR in both the MCA/MSU 
and Korean cattle populations. The 25 lY polymorphism of 
HRG was significantly associated with CW in the Korean 
cattle (p = 0.0167). The IL22-73R polymorphism was 
significantly associated with BF in the MCA/MSU cattle. 
No BTA5 SNP marker was associated with average daily 
gain (ADG) although BL4, BR2936, and BMS1248 MS 
alleles on BTA5 differed in frequencies between high and 
low ADG groups. 

DISCUSSION 

Finding genes responsible for economically important 
traits is difficult because such traits are controlled by 
quantitative genetic factors compounded with 
environmental influence. We hypothesized that major 
genetic factors for economically important traits in beef 
cattle would be detectable when investigating QTL regions 
across cattle breeds. Therefore this study evaluated genetic 
polymorphisms found in commercial cattle from both Korea 
and the USA. There has been no previous study directly 
comparing QTL regions between Hanwoo and other breeds, 
but it is already known that some of the genetic effects 
responsible for highly heritable traits like growth and 
carcass traits are found to be segregating across breeds as 
evidenced from several experimental crosses in beef cattle 



(Abe et al., 2008). 

Bovine chromosomes 1 and 5 were chosen for 
evaluation because QTL regions on these chromosomes 
were detected in an experimental population derived from 
crossing Angus and Brahman (Kim et al., 2003), and 
segregation of these QTL has been validated in other cattle 
populations (Morsci et al., 2006; Li et al., 2004). In addition, 
microsatellite markers on BTA 1 have been associated with 
carcass traits in Hanwoo cattle (Choi et al., 2006). In this 
study, eight microsatellite markers spanning QTL regions 
on BTA 1 and 5 were genotyped in beef cattle from 
Michigan in the USA. Individual allelic frequencies were 
compared between high and low phenotypic groups and chi- 
square tests were performed to determine significant 
differences between the two groups. These initial results 
from the BTA 1 and 5 QTL regions were encouraging to 
suggest that these chromosomal regions might be associated 
with quantitative phenotypic variation for growth and 
carcass traits in commercial beef cattle populations. 
However, further analysis of the allelic association was 
problematic using microsatellite markers because some of 
the alleles had very low frequencies so sample size for an 
association study was reduced. 

To develop genetic makers that were polymorphic in 
both populations, we screened candidate genes located in 
the QTL regions for polymorphisms using available cattle 
genomic data. A total of 25 fragments from 24 candidate 
genes were sequenced in both Korean and MCA/MSU 
cattle samples. Twenty-seven SNPs were identified and 
Unkage disequilibrium was determined among the 
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polymorphisms. Frequencies of these polymorphic markers 
were compared between the two populations and 
phe no type -genotype association analyses were conducted. 
Of the 17 SNPs for which association analyses were 
performed, 10 SNPs exhibited a significant association with 
at least one phenotypic trait. It was not feasible to construct 
haplotypes between the markers because physical distances 
of these markers were too large and the number of animals 
was relatively small for haplotype construction. Some of the 
associations were detected in both populations, although the 
directions of the genotypic effects were not consistent, 
which may indicate different linkage phase of the maker 
allele with the given trait of interest across populations. 
Also, due to different commercial sources for the 
phenotypic data in the two countries, the traits recorded 
were not consistent and even in cases where similar traits 
were determined, the measurements may have differed. 
Thus, this underscores the difficulty of comparing data 
across populations. However, the results support the 
observation that QTL for growth and carcass traits are 
segregating on BTAl and 5 in both populations. 

Currently high-density SNP genotyping platforms are 
available for several domestic species and genome-wide 
association studies have been completed in several cattle 
populations (Fan et al., 2010). In the future, for populations 
with similar genetic characteristics, QTL identification will 
utilize high-density SNP arrays. It is important to classify 
phenotypes compatibly across breeds, and validation of the 
identified QTL regions using smaller density SNP arrays 
may be cost-effective to use for local livestock breeds such 
as Korean cattle (Han woo). Our results in this study may 
contribute to the development of efficient approaches for 
evaluating DNA maker effects and their utilization for 
marker-assisted selection in local conmiercial breeds. 

IMPLICATION 

In summary, markers in BTA 1 and 5 chromosomal 
regions were genotyped in two conmiercial beef 
populations. These chromosomal regions were previously 
associated with growth and carcass traits in several cattle 
crosses. Significant allelic associations were observed with 
growth and carcass trait measurements in this study, which 
indicates that the QTL on these chromosomes may be 
segregating in multiple conmiercial beef populations. Based 
on the QTL locations, 24 genes were considered as 
positional candidate genes within the BTA 1 QTL region 
between 51.0 cM to 88.42 cM, and the BTA 5 QTL region 
between 61.85 cM to 99.16 cM. Comparative re-sequencing 
identified 27 SNP markers of which 17 SNP were evaluated 
for phenotypic associations in Korean and Michigan (USA) 
commercial beef cattle populations. More work is necessary 
with more genes and animals to utilize the DNA marker 



information on BTAl and BTA5 for improved beef 
production, but the identified and verified SNPs from this 
study provide useful information for development of DNA 
marker panels for commercial application. 
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